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For a class of idealized chaotic systems (hyperbolic systems) correlations decay exponentially in 
time. This result is asymptotic and rigorous. The decay rate is related to the Ruelle-Pollicott 
resonances. Nearly all chaotic model systems, that are studied by physicists, are not hyperbolic. 
For many such systems it is known that exponential decay takes place for a long time. It may 
not be asymptotic, but it may persist for a very long time, longer than any time of experimental 
relevance. In this review a heuristic method for calculation of this exponential decay of correlations 
in time is presented. It can be applied to model systems, where there are no rigorous results 
concerning this exponential decay. It was tested for several realistic systems (kicked rotor and 
kicked top) in addition to idealized systems (baker map and perturbed cat map). The method 
consists of truncation of the evolution operator (Frobenius-Perron operator), and performing all 
calculations with the resulting finite dimensional matrix. This finite dimensional approximation can 
be considered as coarse graining, and is equivalent to the effect of noise. The exponential decay rate 
of the chaotic system is obtained when the dimensionality of the approximate evolution operator 
is taken to infinity, resulting in infinitely fine resolution, that is equivalent to vanishing noise. The 
corresponding Ruelle-Pollicott resonances can be calculated for many systems that are beyond the 
validity of the Ruelle-Pollicott theorem. 



I. INTRODUCTION 



The purpose of the lectures that are summarized in this review is to describe the behavior of ensembles of chaotic 
systems. As in the case of statistical mechanics, the dynamics of ensembles turns out to be simpler then the one 
of individual systems. For chaotic systems Jl]-[7j the long time asymptotic behavior is rigorously known to exhibit 
exponential decay of correlations only for a class of idealized systems (hyperbolic systems) . Physical model systems do 
not belong to this class. In what follows nonrigorous methods, that enable the exploration of the long time behavior 
of chaotic systems, are presented and their application is demonstrated for several systems. The review is pedagogical 
and descriptive in nature, and is intended for an overview of the subject. The reader should consult the references 
for the details and the precise statements. 

The dynamics of the systems considered in this review is of Hamiltonian nature. The dynamics of the continuous 
systems is determined by the Hamilton equations: 



V - dq 

a= 22* 
i dp 



(1) 



where Ti is the Hamiltonian, while q and p are the position and momentum respectively. The phase space density 
p(x), where x = (q,p) are the phase space points, satisfies the Liouville theorem, 

f = ° < 2 > 

or 

%={^P}^ (3) 

where {...} are the Poisson brackets. The main issue that will be discussed is the way p spreads in phase space for 
chaotic systems. For classical chaotic systems finer and finer structures develop for longer and longer times. These 
structures are reflected in p. Coarse graining, with some fixed scale of resolution, results in the truncation of the 
evolution of these very fine structures. On the finite (but arbitrarily small) resolution scale the equilibrium uniform 
phase space density is approached. In this review the asymptotic (in time) relaxation to this density will be discussed 
and in the end the coarse graining scale will be taken to 0. The results differ from the ones obtained without coarse 
graining, since the limits of infinite time and vanishing coarse graining do not commute. The approach, where a 
finite resolution scale is used and then the limit where this scale tends to zero is taken, is relevant for experimental 
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realizations where the idealized classical description on the finest scales is destroyed resulting of the coupling to the 
environment. 

Maps are transformations of phase space in discrete time, denoted by n. A map F is a transformation 

x„+i = F(x„). (4) 

We study maps since they are easier to handle analytically and numerically and they reproduce the most interesting 
results found for systems evolving continuously in time. 
Maps can be derived from Hamiltonians of the form: 

H = ^ + V(q)J26(t-n). (5) 

n 

For such maps the phase space area is preserved. Maps can also be defined with no reference to a Hamiltonian. The 
dynamics of phase space densities will be explored here for several area preserving maps: a. Kicked Rotor (Standard 
Map), b. Kicked Top, c. Arnold Cat Map and d. Baker Map. 

a. Kicked Rotor (Standard Map) 
A planar rotor that is periodically kicked is modeled by the Hamiltonian 

72 

H = — + K cos 9 J25{t-n), (6) 

n 

where 9 is the coordinate and J is the conjugate momentum. The Hamilton equations ([!]) are: 

9 = J (7) 
j = K sin 6 ^2 S(t-n). 

n 

Integration with respect to time results in the Standard Map 

@n+i = &n + Jn (8) 
Jn+l = Jn + Ksm6 n+ i, 

where 9 n and J n are the angle and angular momentum just after the n-th kick. The equations (0) define a map of 
the form (||) 

{9 n +i, J n +i) — F(6*„, J n ). (9) 

It is easily checked that the map is area preserving. It becomes more chaotic as the stochasticity parameter K 
increases. The phase space is plotted in Fig. |l| 

A variant of this map provides a good description of driven laser cooled atoms § and of beams deflected by 
dielectrics with a modulated index of refraction HI . 
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FIG. 1. Phase space portraits of the weakly chaotic kicked rotor with stochasticity parameter K = 2 (a), and of the strongly 
chaotic kicked rotor with K = 10 (b). 



b. Kicked Top 

A large spin can be described by a classical vector 

J = j (sin 9 cos ip, sin 9 sin <p, cos 9) 
where 9 and ip are the polar angles. The kicked top map is defined by the transformation (|) 

F = R z {Tcos9)R z ((3 z )R y ((3y), 



(10) 



(11) 



applied to J where Ri(P) is the rotation around the axis i by the angle (3. The nonlinearity results of the dependence 
of a rotation on the angle 9. The chaoticity of the map increases with r that is the stochasticity parameter. The 
canonical phase space variables for this map are q = ip and p = cos 9. The phase space is plotted in Fig. 0. 
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FIG. 2. Phase space portraits of the weakly chaotic kicked top with stochasticity parameter r = 2.1 (a), and of the strongly 
chaotic kicked top with r = 10.2 (b). (Fig. 1 of @). 



c. Arnold Cat Map 

It is defined in the [0,1] x [0, 1] square of the (or, y) phase plane by j^fl] 

%n+l = Xn + Un mod 1 

y n+1 = x n + 2y n mod 1. 

The evolution is demonstrated in ^,0 . 

d. Baker map 

It is defined in the [0,1] x [0, 1] square of the (x, y) phase plane by jijjfjjl] 



(x n+ i,y n+ i) = F(x n ,y n ) = 
This transformation is demonstrated in Fig. |^. 



(2x„,y n /2) 

(2x n - 1, (y n + l)/2) 



(12) 



for < x < \ 
for \ < x < 1 



(13) 
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FIG. 3. The baker map. 



The systems (a) and (b) model physical problems. Their phase space is mixed, namely the dynamics in some parts 
is chaotic and in some parts it is regular. Systems (c) and (d) are very idealized. For these systems the motion is 
chaotic in the entire phase space. 

A system is chaotic if points, that are initially close, spread exponentially in phase space. For chaotic systems the 
separation between points grows with the number of steps as 

e(n) = e(0)A", (14) 

in the limit n — > oo and e(0) — » 0. The Lyapunov exponent is 

A = In A (15) 

and A is sometimes called the Lyapunov number. A system that evolves in a bounded region in phase space is called 
chaotic if its Lyapunov exponent A according to the definition ( ^|Jl5| ) is positive. (More generally the Lyapunov 
numbers are the eigenvalues of the tangent map. The Lyapunov number defined by (|14| ) is equal to the magnitude of 
the largest eigenvalue.) One can study the local expansion and contraction in the various directions. If for a map, for 
any point of phase space, in each direction, either expansion or contraction takes place, or in other words if none of 
the real parts of the local Lyapunov exponents vanishes, the system is called hyperbolic. 

Exponential spread of trajectories in phase space does not imply exponential decay of correlations in time. Such 
correlations may result of sticking to small structures in phase space, as is the case for the stadium billiard and for 
the Sinai billiard as well as for chaotic trajectories of mixed systems, such as the kicked rotor and the kicked top. For 
hyperbolic systems, such as the baker map, typically correlations in time decay exponentially. 

Chaotic motion is reminiscent of a random walk. Therefore it is instructive to review the dynamics of the probability 
density of random walkers. Assume a random walk on a one dimensional lattice, with equal probability of 1/2 to 
make a step to the left or to the right. The path of each walker is complicated but the evolution of a probability 
density of random walkers is simple. It is simple in particular in the continuum approximation (the limit of vanishing 
lattice spacing) where it is described by the diffusion equation 

-p(x,t) = D-^p{x,t), (16) 

D is the diffusion coefficient and p{x, t) is the density of random walkers at the point x at time t. Let us assume that 
the random walk is confined to an interval of length s, for example —s/2 < x < s/2. The current at the ends of the 
interval vanishes, resulting in the boundary conditions 



dp 

dx 



X = ±s/2 = (17) 
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on ([l6]) . This problem is similar to the standard quantum mechanical problem of a particle in an infinite square well, 
but with the boundary conditions ( p"7| ) and with no ih on the left hand side of ([if]). Therefore the density of states 
can be expanded as: 



where 



Uk 



form an orthonormal basis and 



p(x,t) =^a fc e lkt u k (x), 



fc=0 



sin — k 



2 

s s 



(18) 



k even 
k odd 



(19) 



(20) 



The units of a k are of 1/ ^/[length]. The expansion coefficients a k are determined by the initial density p(x,t = 0). 
The equilibrium density is 



lim p(x,t) = ao\ -, 

t^co V S 



and the asymptotic approach to this density is 

p(x,t) 



a\e 7lt Ui(x) 



(21) 



(22) 



with the relaxation rate 7i = (7) D. The density-density correlation function is 

rs/2 



dxp(x, 0)p{x, t) - a 2 Q = V 

-s/2 



(23) 



fc=i 



Hence also the correlations decay with the rate 71 . 

The example of the probability density of random walkers demonstrates that although the trajectories of specific 
random walkers are complicated their probability density p follows simple dynamics. We turn now to study the phase 
space probability density p(x), that for Hamiltonian dynamics satisfies (||). The introduction of such probability 
densities, that are sufficiently smooth, is actually a coarse graining over some fine scale. It is an averaging over phase 
space of the same type that was performed for random walkers. For any given area preserving map F it is instructive 
to introduce the one step evolution operator U , 



p„+i(x) = Up n (x), 



(24) 



where /0«(x) is the phase space density after n steps of the map. If the operator is defined on a space of sufficiently 
smooth functions it is called the Frobenius-Perron operator. If F is area preserving and invertible, namely if F _1 
exists, then U is unitary. All eigenvalues of U are on the unit circle in the complex plane, they take the form e ^ i27rQ j 
with real a. The corresponding "eigenfunctions" are not square integrable. An example of such a function will be 
presented now mW. Assume the map F has a periodic orbit of period n, that is: 



x i+i 
xi 



F(x,), 
F(x„). 



3 



1,2, ...n 



(25) 



The function 



tA q (x) =]T ( 5(x-x J )e^ 



(26) 
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with a — l/n, where I is an integer, is an "eigenfunction" of U with the eigenvalue e l2va . Here ip a is called 
"eigenfunction" if 

Ui> a = e- i27 "*ip a (27) 

and no requirement is made that it belongs to the space where U is defined. For example if U is defined on the space 
of square integrable functions, ip a of (|2^) does not belong to this space and is not a function in the usual sense (it is 
a distribution). 

The operator U can be represented by an infinite dimensional matrix. In physics applications it is natural to use 
finite dimensional approximations. Let us assume that a basis that is ordered by increased resolution is used (the 
basis states may be for example sines and cosines or orthogonal polynomials). For the truncation of U , namely a finite 
dimensional approximation, the matrix is not unitary, and its eigenvalues are inside the unit circle in the complex 
plane. They may vary with the dimension of the matrix. The natural question is: Do these eigenvalues approach 
the unit circle in the limit of an infinite dimensional matrix? It will be demonstrated in what follows that this is not 
the case for chaotic systems and values inside the unit circle are approached in this limit, in spite of the fact that 
the infinite dimensional matrix representing U is unitary. A heuristic justification for this behavior was proposed by 
F. Haakc Chaotic systems, in contrast with regular ones, exhibit phase space structures on all scales. These 

structures are revealed during the evolution. The operator U couples the fine scales via its matrix elements that 
couple to states with high resolution. As a result of the truncation, probability that was originally transferred to 
the fine scales is lost, resulting in nonunitarity of U . Convergence of the eigenvalues to values inside the unit circle 
in the complex plane, in the limit of infinite dimension of the matrix, results of the asymptotic self similarity of the 
chaotic dynamics. For regular motion, on the other hand, as the dimension of the matrix is increased the eigenvalues 
approach the unit circle (see Sec. IV). 

It will be shown that for chaotic systems the finite dimensional approximations of U, in the limit of large dimension, 
describe the decay of correlations. The phase space density-density correlation function, in analogy with (E3[) is 



C(n) = J dxp(x, 0) P (x, ti) - p^n (28) 

= J dxp(x,0){J>(x,0)-p^fi, 

where fi is the volume of the chaotic component in phase space. The equilibrium density is = lin^^oo p(x, n) 
that is independent of position in phase space. A more general correlation function is: 

C^ B \n) = J dxA(x)L>" J B(x), (29) 

where we assumed 

lim U n A{x) = lim U n B(x) = 0. (30) 

n — >oo n — >oc 

This can always be obtained by the subtraction of the asymptotic value. For simplicity it will be assumed that both 
A and B are real. 

It is useful to study the Laplace transforms of the correlation functions. For this purpose we introduce the resolvent 

oo 

R(z) = - V U j z~ j = (31) 

z ^ z — U 

3=0 Z U 

where z is a complex number. Since U is unitary the sum is convergent for \z\ > 1. This is analogous to the usual 
definition used in quantum mechanics: 

1 f°° - (% e 

<>■"■■■ = " 

where 



R(E) = — / dtUexp -Et --rt], (32) 
in J \n n 1 



U = exp[-^Ht) (33) 
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is the evolution operator, leading to 



R(E) = 1 ■ (34) 
E-H + it 



for e > 0. The convergence of the integral (|32|) requires e > 0. In analogy the convergence of the sum (31) requires 
\z\ > 1 in (EH]). The discrete Laplace (one sided Fourier) transform of the correlation function (p9t) is 



& A > B \z) = jr C {A > B \n)z- n = [ dxA(x) \zR(z)} B(x), (35) 

n=0 ^ 



as one finds from (pl|). 

Let V'i be an "eigenfunction" of U with the eigenvalue Zi 



Uipi = Ziipi. (36) 



By "eigenfunction" we mean here that it satisfies (pq), but it may not be a function in the usual sense (for example 
it may be a distribution) . Then Zi is a pole of the matrix elements of the resolvent as one can see from (|3l]) . The 
correlation function, involving an "eigenfunction" , takes the form 

C^(z) = (A\i>i)^—, (37) 

Z - Zi 

where the Dirac notation 

(A\B) = I dx J 4*(x)B(x) (38) 



is used (A* is the complex conjugate of A). Here A and B were assumed to be real. Let us expand B(x) in terms of 
the "eigenfunctions" of U as 

B(x)=5^6^i(x). (39) 

i 

Then 

C^ B \z) = -D— (40) 
z- z l 

where the poles Zi depend only on U while the residues 

n = bi{A\i/ji)z, (41) 

depend also on A and B. 

Because of the unusual nature of the "eigenfunctions" tpi the manipulations leading to ([hi]) are only heuristic. The 
Ruelle-Pollicott theorem justifies the expression ( fib] ) for hyperbolic systems if A and B are typical smooth functions 
(for a precise statement of the theorem that is transparent for physicists see [0). Moreover the theorem assures that 
\zi\ < 1, except for the eigenvalue zq = 1 corresponding to the equilibrium density ipo, that is independent of x. The 
Zi are called the Ruelle-Pollicott resonances. 

The existence of poles with \zi\ < 1 implies the decay of correlations with the rate In This results of 

C^( n ) = — I ^lc^ B \z)z n (42) 

and the application of the residue theorem to (|40|). 

For any approximation of dimension TV of the Frobenius-Perron operator U, that will be denoted by U( N > in what 
follows, the eigenfunctions i>\ , corresponding to the eigenvalues z\ N \ are well defined. For hyperbolic systems, 
where the Ruelle-Pollicott theorem applies, it is reasonable to assume that in the limit N — > oo, the eigenvalues z\ N ^ 
approach the Ruelle-Pollicott resonances, for a typical choice of the basis. Inspired by the Ruelle-Pollicott theorem 
the following heuristic scheme for the calculation of the decay rates of correlations is proposed: 
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1. Introduce an orthogonal basis where the basis states are ordered by resolution. These states, for example, may 
be orthogonal polynomials or trigonometric functions (cos lukx, sin 2nkx) or exponentials exp i2irkx. 

2. Calculate the matrix elements of U in this basis. 

3. Introduce a truncation of dimension N of this matrix, that will be denoted as U^ N K 

4. Calculate the eigenvalues z\ N ^ and the eigenfunctions ip\ N ^ of U^ N \ The eigenf unctions ^j- are finite linear 
combinations of the basis states, and therefore are smooth. 

5. Take the limit N — > oo. In this limit zj- — * Zi and ip\ N ^ — > Vi- 

From the Ruelle-Pollicott theorem it is expected, for hyperbolic systems, that the Ruelle-Pollicott resonances are 
obtained in this way, if a typical basis is used. The limiting functions ipi are singular |13|-|l5[j . This scheme was 
applied to various systems, including also to mixed systems |jl(j , [Lg -g(j. Limiting eigenvalues z\ satisfying \zi\ < 1 
were obtained in this way, implying exponential relaxation of correlations with the rate ln|zj| that takes place for a 
very long time (for some nonhyperbolic systems it is known that eventually the correlations decay as a power-law, 
but the time required to obtain this power-law behavior may be too long for any physical relevance). Although the 
Hamiltonian dynamics of trajectories is reversible the phase space densities exhibit relaxation that is irreversible, and 
is similar to the behavior of probability densities of random walkers that follow irreversible dynamics. 

In Sec. II the heuristic scheme for the calculation of the Zi will be demonstrated for the baker map that is hyperbolic 
and consequently the Ruelle-Pollicott theorem applies. Exploration of a hyperbolic system that is a modification of 
the cat map is briefly mentioned in Sec. III. In Sec. IV the scheme will be applied to the kicked top and in Sec. 
V it will be applied to the kicked rotor, that are mixed systems, where the Ruelle-Pollicott theorem is not valid. 
Also the limitations on the validity of this scheme for such systems will be demonstrated. The main conclusions 
are summarized in Sec. VI. The review will follow closely references 1 10 L^|nj, that are marked by ** in the list of 



references. A detailed discussion on the relaxation of correlations in chaotic systems and on related topics can be 
found in • The scheme for the use of truncation in the calculation of the Ruelle-Pollicott resonances is discussed 
in a wider context in pTj . 

II. THE FROBENIUS-PERRON OPERATOR FOR THE BAKER MAP, A DEMONSTRATION 

The heuristic prescription for the analysis of the Frobenius-Perron operator that was outlined in the Introduction 
will be demonstrated for the baker map (|l3|), where all the results are exactly known. The analysis will make use of 
exact results obtained in 

The basis states that will be used are 

(xy\kl) = P k (x)P l (y) (43) 

where 

P k (x) = V2k + lPk(l - 2x) = V J ^ x k (l~x) k (44) 

are the modified Legendre polynomials while P k are the Legendre polynomials. A basis of orthogonal polynomials is 
natural, since the order of the polynomial is not affected by the map (|l3|). The orthonormality of the basis 

(kl\k'l') = 5 kk >8 w (45) 

follows from the orthonormality of the modified Legendre polynomials 

l 

dxP k {x)P k ,{x)=5 kk ,. (46) 

The basis is naturally ordered by increased resolution, since this is a property of orthogonal polynomials. The action 
of U of ( p4] ) on any phase space density is 

Up(x) =p(F- 1 (x)), (47) 
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since the points that are at x at time step n + 1 were at F 1 (x) the time step n. For the baker map ( |l3| ) 

p(F (j,(( x+1 y 2i 2y-l) for±<y<l. (48) 

From @ and @ one finds that the matrix elements of U are 

(kl\U\k'l') =\[l+ (-l) fc+fc '+ w '] J^,/,,, (49) 

where 

I kk , = / dxP k {x)P k ,{x/2). (50) 
Jo 

Using (Q) and integrating (f50| ) by parts k times one finds that if k > k', 

hk> = 0. (51) 

This results in the nonrecurrence property of U 

(kl\U\k'l') = for k > k' or I' > I. (52) 

During the evolution, probability is transformed from states with k! to states with k only if k < k' , therefore 
after application of U the density becomes more uniform because the weight of lower order Legendre polynomials is 
increased. This is expected since x is the unstable direction, where stretching takes place, making the density more 
uniform. In the y direction, on the other hand, the density is transformed from /' to I only if I > therefore during 
the evolution, the weight of the high order Legendre polynomials increases. This results of the fact that y is the stable 
direction, where contraction takes place, resulting in complexity that increases with time. 
For k < k' one finds |TJ 

[(2fc + l)(2fc- + l)] 1 / 2 ^7 lV (fc' + fc + Q! 

kk ' 2 k t^V 2/ (k'-k-l)\(2k + l + l)HV [ ' 

In particular I kk = 2~ fc and the diagonal matrix elements are 

(kl\U\kl) = JL_. (54) 

The matrix elements of U n can be calculated with the help of the nonrecurrence property (^). Introducing the 
resolution of the identity one finds 

(kl\U n \k'l') = (kllUlhl^ihhlUfah) (kn-lln-l\UWl'), (55) 

[**;/<] 

where [fcj;Zj] = [fci, k2 fcn-ij ^1; fa, ■•■■ln-i\ an d the nonrecurrence property implies k < k\ < &2 < k 1 and 

I > li > I2 > V ■ In particular the diagonal matrix elements satisfy 

(kl\U n \kl) = (H\U\kir = (j!^) ■ (56) 

In the limit n — > 00 the off-diagonal matrix elements are dominated by powers of the diagonal matrix elements, since 
the nonrecurrence property limits the number of non-diagonal matrix elements in (|55"|). 

The diagonal matrix elements of the resolvent can be easily calculated with the help of (|3l|). For \z\ > 1 the sum 
in (31) is convergent and using ( |56| ) one finds 

00 ^ 
(kl\R(z)\kl) = £ (kl\U\kl) n = z _ 2 _ (k+iy (57) 

n=0 
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These matrix elements are singular for z — 2~( fe+ ') , and except for k — I — all these singular points are inside the unit 
circle in the complex plane. This is also the case for the off-diagonal matrix elements. Therefore the Ruelle-Pollicott 
resonances are 

z m = 2~ m (58) 

and their degeneracy is m + 1. They are related to the decay ol correlations by ([lO]) . They were obtained by analytic 
continuation from \z\ > 1 where the resolvent is defined. One can show that the matrix elements exhibit a cut at 
\z\ = 1 and the physically relevant poles are on the Riemann sheet continued from \z\ > 1. 

What happens if the Frobenius-Perron operator is restricted to k, k' < k max and I, I' < l m ax and U is approximated 
by IJ( N \ an N = k max l max dimensional matrix, resulting of the truncation of Ul The nonrecurrence property holds 
for L/W and consequently ( p6| ) holds. Therefore for the N dimensional matrix f/W the diagonal matrix elements, 
are zeros of the characteristic polynomial, that is of order N. To see this, note that by the Hamilton-Cayley theorem, 
satisfies its characteristic polynomial, namely ^^-qPjU^^ = 0, and calculate the diagonal matrix elements of 

this expression with the help of ([56]) . Consequently the diagonal matrix elements are eigenvalues of U( N ^ , taking the 
values 2~ m . The multiplicity of the eigenvalue 2~ m is m + 1 if m < k max , l max . By a similarity transformation the 
matrix can be transformed to the canonical Jordan form. Another way to see that the eigenvalues of U^ N ^ coincide 
with the diagonal matrix elements is by using the indices l max — I and l max — I' instead of I and I'. With these indices, 
is upper triangular, therefore the eigenvalues coincide with the diagonal matrix elements. The eigenvalues are 
independent of k max and l m ax, therefore these are of the form 2 _m also in the limit k max — > oo and l max — > oo, in 
spite of the fact that U is unitary. For U( N > the eigenfunctions are finite combinations of Legendre polynomials and 
are therefore smooth. In the limit k max — > oo and l max ~* oo the right eigenfunctions are independent of x and are 
polynomials of y, while the left eigenfunctions tend to distributions jOl. 
The main results that were obtained for the baker map are: 

1. A basis that is ordered by increased resolution was introduced. A truncation of dimension in this basis was 
implemented. 

2. The Ruelle-Pollicott resonances were calculated from U^ N \ the N dimensional truncation of U . As N increases 
more eigenvalues are revealed. The values of the eigenvalues do not depend on the truncation dimension A, and 
therefore also in the limit of infinite N they remain at values that were obtained for finite N . 

3. The poles of the matrix elements of the resolvent (|3l| ) and ( |57|) were obtained by analytic continuation from 
|*| >1. 

4. The eigenfunctions of U are functions of y only (y is the stable direction). 

In this section the procedure for the calculation of Ruelle-Pollicott resonances from JJ^ N \ a truncated matrix 
approximating U, was demonstrated for a system where the results are known exactly. In the following sections this 
method will be applied to systems where there is no exact theory. 

The basis used here was of Legendre polynomials. If a basis of sines and cosines is used it is impossible to find the 
Ruelle-Pollicott resonances in this way, since each of the basis states collapses to after a finite number of applications 
of U |l^j2l[] . On the basis of the Ruelle-Pollicott theorem, we believe that for a typical basis one does not encounter 
problems of this nature. 



III. THE FROBENIUS-PERRON OPERATOR OF A MODIFIED CAT MAP 

For the Arnold cat map, defined by (|l2|), correlations in time decay faster than exponentially. A modification, 
where the function /(x„ + i), that is defined by 

f(x) = — - [cos(27re) — cos(47rx)] , (59) 
2n 

is added to the equation for y n +i in (|lj), is a hyperbolic system if Kq is sufficiently small [Q. The correlations for 
this system decay exponentially in time. It was studied with the help of the truncated Frobenius-Perron operator, 
according to the scheme that was outlined in the introduction. The leading Ruelle-Pollicott resonances were found 
with the help of a variational approach. 
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IV. THE FROBENIUS-PERRON OPERATOR FOR THE KICKED TOP 



We turn now to apply the method where U is approximated by a finite dimensional matrix to mixed systems (where 
in some parts of phase space the motion is chaotic and in other parts it is regular). The Ruelle-Pollicott theorem does 
not apply to such systems. It will be demonstrated, however, that the Ruelle-Pollicott resonances are meaningful 
and describe the decay of correlations in the chaotic component of mixed systems for a time that may be very long. 
Asymptotically power-law decay takes place, due to sticking to regular islands in phase space. 

In this section the kicked top H, that was defined by (10) and (fy]), will be analyzed. It is a summary of results 
that were presented in |lO| (see also |l7j] and p8[). The phase space is shown in Fig. |2[ 

A natural basis for this problem is of spherical harmonics 



Y lm (9, (p ) = (-iy 



'21 + 1 (l-m)\ 
47T (l + m)\ 



P ; m (cos( 



(60) 



Alternatively one can choose a basis of real functions with e lmv replaced by smrrup and cos my, resulting in a real 
matrix for U. It is naturally ordered by increased resolution. A truncation at I = l ma x is introduced. The dimension 
of the truncated matrix U^ N \ approximating U is N = (l max + l) 2 . As l max is increased the eigenvalues of U^ N ' 
converge to values inside the unit circle in the complex plane, as can be seen from Table Q. 



^max — 30 


/max = 40 


^max — 50 


^rnax — 60 


0.7700 


0.7688 


0.7523 


0.7696 


0.3075 


0.3429 


0.3523 


0.3550 


±i 0.5740 


±i 0.6140 


±i 0.6211 


±i 0.6199 


-0.3170 


-0.3348 


-0.3444 


-0.3388 


±i 0.6003 


±i 0.6272 


±i 0.6283 


±i 0.6243 


-0.0042 


-0.0002 


-0.0100 


-0.0058 


±i 0.7161 


±i 0.7133 


±i 0.6930 


±i 0.7080 


-0.7025 


-0.7228 


-0.7155 


-0.7165 


0.6544 


0.6230 


0.6495 


0.6480 




-0.5619 


-0.5753 


-0.5667 



TABLE I. Eigenvalues of U (N) for r = 10.2 truncated at Z max = 30, 40, 50 and 60 (Table I of Eol). 
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The corresponding eigenfunctions are presented in Fig. ^. 




(b) 



-1 




-11 




FIG. 4. The eigenfunctions corresponding to the eigenvalues 0.7696 (a), -0.3388±i0.6243 (b), -0.0058±i0.7080 (c), 0.6480 
(d) of U (N) for r = 10.2, and Z max = 60 (Fig. 4 of @). 

Dark-shaded regions in phase space indicate large amplitudes of the eigenfunctions. In Figs. || (a) and (b) 



(a) 





(c) 



m 



-IT 



q it 

FIG. 5. As phase-space resolution increases from Z max = 30 (a) to Z max = 60 (b), the eigenfunction of for the eigenvalue 
—0.3388 ± i0.6243 (for r = 10.2) gains new structures on finer scales. The corresponding eigenfunction (c) of [/ T ( iv ' (resolution 
'max = 60) is localized at the same periodic orbit as the eigenfunction (b) but with stable and unstable manifolds interchanged. 
See also figure |^(a) for the periodic orbits, figure ^|(b) for the unstable and figure ^|(c) for the stable manifolds (Fig. 5 of jnj). 
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it is demonstrated that as l max increases finer details of the eigenfunctions are revealed. This is a result of increased 
resolution. Comparison with Fig. [| (b) 





FIG. 6. (a) : The 12 orbits of primitive length 6 that can be identified for the resonances 0.3550±i0.6199 and -0.3388±i0.6243 
related to the eigenfunction shown in figures^(b) and ^(b). The unstable manifolds of these orbits are shown in (b), and the 
stable manifolds in (c) (Fig. 6 of Jlo|). 

demonstrates that the eigenfunctions are large and tend to be uniform on the unstable manifold. In Fig. [?] 

C(n) 

0.6 
0.4 




5 10 15 n 

FIG. 7. The decay of C(n) (dots) with p(0) corresponding to the eigenfunction shown in figure H(b) (also plb)). The 
numerical fit (line) yields a decay factor 0.7706 (compared to l-z'^'l ~ 0.7103). (Fig. 7 of Jlo|). 

the density-density correlation function (proportional to (p8|)) is plotted. It is demonstrated that also here the 



14 



Ruelle-Pollicott resonances describe the decay of correlations for a very long time. 

The behavior of the eigenfunctions localized in regular regions and the corresponding eigenvalues is very different. 
Eigenfunctions that are localized in islands of regular motion presented in Fig. || (a) are depicted in Fig. ||. 




P 



FIG. 8. Eigenfunctions of with Z max = 60 that are localized on elliptic islands shown in Fig. 2 (a) for the weakly 

chaotic case r = 2.1. The corresponding eigenvalues 0.999976 (a) and 0.999974 (b) are almost at unity (Fig. 3 of fic|l). 



The eigenfunctions follow elliptic orbits inside the regular island and therefore do not involve the resolution of very 
fine structures. The corresponding eigenvalues approach the unit circle in the complex plane as l m ax increases. 
It is instructive to examine the evolution operator U^ 1 corresponding to the inverse map F _1 . Unitarity implies 

U~ 1 = Ul (61) 

A basis can be chosen so that U is real and W — U T , where T denotes transpose. The eigenvalues of U T are equal 
to the ones of U. Introducing the truncation one hnds 

(jT{N) = £(JV)T (62) 

where {J 7 ^- 1 is the truncated U T and U^ 7 ' is the transpose of U^ N \ the truncated U. The eigenfunctions of U T ( N "> 
are localized on the unstable manifold of F -1 , that is the stable manifold of F, as is clearly demonstrated comparing 
Fig. ||(c) with Fig. || (c). By eigenfunctions we meant so far the right eigenfunctions. The left eigenfunctions of U 
(and U^) are the right eigenfunctions of U T (and U T ( N ^). 
The main results that were obtained for the kicked top are: 

1. As N or l max increases the eigenvalues of U^ N \ corresponding to eigenfunctions localized in the chaotic compo- 
nent of phase space, converge to values inside the unit circle. Unlike the case of the baker map, the eigenvalues 
z\ N ' depend on N and Umjv ->«> zjf^ = Z\ with \z\\ < 1, except the unit eigenvalue corresponding to the 
equilibrium density. 

2. These eigenvalues determine the decay of correlations in the chaotic component of phase space. 

3. In the chaotic component of phase space the right eigenfunctions of the truncated Frobenius-Perron operator 

are localized on the unstable manifold of F, and tend to be uniform along this manifold. The left 
eigenfunctions are localized on the stable manifold. These manifolds are not as simple as for the baker map. 
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4. The eigenvalues of corresponding to eigenfunctions localized in regular regions approach the unit circle as 



V. THE FROBENIUS-PERRON OPERATOR FOR THE KICKED ROTOR 

In this section the decay of correlations in time will be studied for the kicked rotor, defined by the Hamiltonian 
(0), and its dynamics is given by the standard map (@) |l]J|,fi]]. The phase space is presented in Fig. [|. It is a mixed 
system (where in some regions of phase space the motion is chaotic and in other regions it is regular). In this section 
we summarize results that were presented in jl9| (see also [^0| , |2j| and jl(| ) . 

For sufficiently large values of the stochasticity parameter K (that takes a typical value), the spread of angular 
momentum is diffusive, to a good approximation, namely 

{{Jn - Jo) 2 ) « 2Dn, (63) 

for large n. The average is over the initial points and is denoted by ( ). A very crude way to get this result is to 

iterate fl) n times to obtain 

71 

J n - J = K sin 6i, (64) 
t=i 

then squaring and averaging with the assumption of absence of angular correlations 

1 

— ( 

2 



(65) 



one obtain s (p3[ ) with D — \K 2 . A more careful calculation, that takes some of the angular correlations into account 
results in |24[] 

D(K) = ^-(l-2J 2 (K) + ...). (66) 



It is actually an expansion in powers of 1/ v K. Because of the mixed nature of phase space it is expected that diffusion 
may not take place asymptotically in time pq |. Moreover, there are values of K (near integer multiples of 2tt) where 
acceleration in momentum is found for some initial conditions, resulting in ((Jn ~ Jo) 2 ) ~ n2 - I n what follows the 
Frobenius-Perron operator will be used to study whether correlations in angular momentum decay as expected for 
a true diffusion process, described in the Introduction (see (|l6|)-(|23])). We will analyze also the decay of angular 
correlations. 

The Frobenius-Perron operator, corresponding to the map (|^) on the torus 

(0 < J < 2ns) (67) 

(0 < 6 < 2tt), 

where s is integer, will be studied. For the usual kicked rotor s — > oo. Noise is added to the standard map by the 
addition of a random variable £ n with variance a 2 on the right hand side of the equation for 6 n +i in (||) following 
f2~ij| . The noise leads to truncation of the Frobenius-Perron operator. The natural basis for the analysis is 

(J9\km) = j— jL exp(zmfl)exp (i— J . (68) 



where k and m are integers. Note that the functions (J9\k, m = 0) form the basis of eigenstates of the diffusion 
operator in the angular momentum J. In this basis, in presence of noise, the Frobenius-Perron operator is 

(km\U^\k'm') = J m ^ m , {^f\ ex P (~Y m2 ) S k-k',ms- (69) 

For a the operator is not unitary. The effective truncation is m < l/cr, \k — k'\ < s/a, \m — m'\ < k'K/s. The 
(7 — > limit , corresponding to the TV — > oo limit of the previous sections, will be taken in the end of calculation. The 
analysis will be performed for 
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s > 1 (70) 

if>l 
kK 

< 1 

s 



and the limit a — > will be taken in the end of the calculation. 

The Ruelle resonances are identified from the matrix elements of the resolvent (|3l]) . It is useful to define a variant 
of the resolvent 

oo 1 

R'(z) = V U j z j = * (71) 

^ 1 - zI7 



that is related to R(z) via 



)=£(*). (72) 



z \z 



Also the matrix elements of R and R! satisfy (|7^) . Let 

i?i2 = (fcimil-HCOIfema) (73) 

and 

i?' 12 = (fc 1 m 1 |i?'(z)|fc 2 m 2 ) (74) 

be the matrix elements of R and R! respectively. Therefore a singularity of R\ 2 at z c implies a singularity of R 12 at 
l/z c and vice versa. The series (31) for R(z) converge for \z\ > 1 while the series ( |7l| ) for -R'(z) converge for \z\ < 1. 
Therefore the Ruelle-Pollicott resonance (located inside) that is closest to the unit circle corresponds to the singularity 
of R 12 that is closest to the unit circle (located outside). This singularity is the radius of convergence of the series 

oo 

where 

dj = (k 1 m 1 \U j \k2m 2 )- (76) 
According to the Cauchy-Hadamard theorem (see pq] ) the inverse of the radius of convergence is given by 



r 



and asymptotically 

const. 



1 = lim sup { \a,j\. (77) 



(78) 



Since the radius of convergence is the singularity of R' 12 (z) that is the closest to the unit circle it satisfies r = 1/ z c 
and ^fa] — > z c in the limit j — > oo (at least for some subsequence of {j}). 

The limit of the series ([76]) was found in the leading order in perturbation theory in the small parameters implied 
by (f70|). The calculation can be performed separately for the various values of fc, resulting in 

/ k 2 K 2 / 2\\ 
z fc = cxp -— h-2J 2 (K)e-° ) ), (79) 



, 4s 2 . 

corresponding to diffusion modes in presence of noise. Taking the limit a — > one finds 

z k = e-^WIs* (80) 
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where D(K) is given by (|66|). This relaxation is similar to the one that is found for usual diffusion (see ( |16| - |23| )). The 
relaxation rates are 



7fc = %D{K\ 



(81) 



that correspond to (ppf). These describe the decay of correlations in angular momentum J. The equilibrium density 
corresponds to Zq. Exploring the subspace involving \k — 0, m), one finds that the Rucllc-Pollicott resonances, 
corresponding to the slowest mode of the relaxation of correlations in angle, take four values 



± iz 



with 



Taking the limit a — > one finds 



~z= Vk2m.(m*fO|cxp(-<7 2 m* 2 /2). 



V\J2m*(m*K)\, 



(82) 
(83) 
(84) 



where m* is the integer m that maximizes | Ji m {mK)\ for a given value of K. A similar result was found in |lq| by a 
variational approach. The corresponding relaxation rate is 

j = \nz. (85) 

In order to test the analytical results, the correlation function of the form (j2£ 

C fg (n) = (f\U n \g), (86) 

where / and g are taken to be basis states of was computed numerically. The results are plotted in Fig. 0. 



0.5p 

°- 
-0.5 

O 

S-1.5- 
-2 
-2.5 



0.5 

-0.5 
-1 

S-1-5 
| -2 
-2.5 
-3 



(a) 



500 1000 1500 2000 2500 3000 



0.5 

Of 
-0.5 
-1 

g-1-5 

|-2 
-2.5 
-3 
-3.5 

" 4 o 



(b) 



1000 



2000 



3000 



4000 







(c) 





























































































1000 



2000 



3000 



4000 







(d) 


~ Q, 




o - ^ 
cr v 

~~ s 









500 



1000 1500 2000 2500 3000 
n 



FIG. 9. The function C fg (n) for: (a) K = 20 ; (b) K = 30 ; (c) K = 40 ; (d) K = 27, for various functions / and g and for 
various values of s (Fig. 1 of fL9[). 
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Clear exponential decay is found. From the slopes, the relaxation rates jk are obtained for the various modes, k. 
The values of D implied by 74, are calculated from (pl|) and compared to the ones obtained from (p(|) in Fig. M. 



3000 r 




120 

FIG. 10. The diffusion coefficient D for K > 10 as found from plots like the ones presented in Fig. 9. Various symbols 
represent results found for different modes k while the solid line represents the analytical value ( p6| ) (Fig. 2 of p^|). 

Good agreement was found. In Figs. ^| and |l0| large values of the stochasticity parameter, K > 10, were used. In 
this regime the regular regions are extremely small (they are invisible in Fig. |l|), therefore the theory presented in 
this review is expected to work very well, as is indeed found in Fig. |l^. In Fig. [Il] the correlation function ( |S6| ) is 
presented for K < 20. 
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FIG. 11. The function C ff An) for: (a) K = 7 ; (b) K 
various values of s (Fig. 3 of ]l9[). 



8 ; (c) K = 3 ; (d) K — 17, for various functions / and g and for 
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From these results the relaxation rates jk were calculated. The resulting values of the diffusion coefficient obtained 
from (81) are compared to the ones found from (66) in Fig. [l^, for K < 20. Here, in contrast to the regime of large 
K presented in Fig. nO, appreciable deviations are found. 
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FIG. 12. The diffusion coefficient D for K < 20 as found from plots like the ones presented in Fig. 11. Various symbols 
represent results found for different modes k while the solid line represents the analytical value (^) . The dashed line represents 
the approximate error, resulting of the truncation of the perturbation theory expansion. The values of D obtained by direct 
simulation of propagation of trajectories are marked by diamonds (Fig. 4 of fill). 

These are related to sticking to regular structures, that is not taken into account in the theory presented in this 
review. To test the decay of correlations in the angle variable, the correlation function ( |86| ) was computed for |/) and 
\g) that are basis states of the form \k = 0, m), and some of the results are presented in Fig. O. 
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FIG. 13. The function C fg (n) for s = 1 with |/) = |01), \g) = [02) and (a) K = 16.3; (b) K = 19.5; (c) K = 12; (d) K = 16. 
The dashed line represents the best fit to the data (Fig. 5 of fl9f| ). 
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The relaxation is very fast, therefore it is very difficult to estimate numerically the relaxation rate. The numerical 
results are compared to the analytical prediction (|5|) in Fig. [li]. 

1.6 
1.4 
1.2 

1 
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0.6 ; 
0.4 
0.2 



"10 12 14 16 18 20 

K 

FIG. 14. The fast relaxation rates y a s found from plots like Fig. [l^ for various functions / and g with k = 0, marked by 
various symbols. The analytical value ( J85| ) is represented by a solid line while the dashed lines denote the analytically estimated 
error, resulting of the truncation of the perturbation theory expansion (Fig. 7 of Jiit). 

The analytical formula provides a reasonable estimate of the relaxation rate and of its dependence on the stochastic- 
ity parameter K. There are however deviations that are significant. Also here some of the most significant deviations 
were found to be related to regular structures. 

The main results that were obtained for the kicked rotor are: 

1. Addition of noise results in the effective truncation of the Frobenius-Perron operator. In the limit of vanishing 
noise the Ruelle-Pollicott resonances are found. The resonances depend on the noise, and therefore on the 
truncation, as is the case for the kicked top, and approach limiting values inside the unit circle in the limit of 
vanishing noise, except the unit eigenvalue corresponding to the equilibrium density. 

2. The Ruelle-Pollicott resonances determine the decay of correlations. 

3. Deviations from exponential decay of correlations may be found after some time as a result of the existence of 
regular structures. 

4. The following physical picture emerges. First the correlations in the angle variable 9 decay. Then the correlations 
in angular momentum J decay as for a true diffusion process. For finite s a uniform equilibrium density is reached. 
For s — oo, after the decay of the angular correlations, diffusive spreading takes place. 

5. At some time, that is increasing with the stochasticity parameter sticking to regular structures becomes 
important, and the picture based on Ruelle-Pollicott resonances and exponential decay of correlations, breaks 
down. 

VI. SUMMARY 

It was demonstrated that Ruelle-Pollicott resonances are relevant even when the conditions for the Ruelle-Pollicott 
theorem do not hold. A method for the calculation of these resonances by a finite truncation of the Frobenius-Perron 
operator was presented. The relation between the resonances and the decay of correlations was established under 
these conditions. For mixed systems (as well as for other nonhyperbolic systems) the picture breaks down after some 
time, that may be very long, as a result of sticking to regular structures. 

In order to evaluate the relevance of the Ruelle-Pollicott resonance picture for realistic systems it is instructive to 
introduce several time scales: 
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1. If 71 is the slowest decay rate (corresponding to the resonance that is closest to the unit circle), t c haos = I/71 
is the decay time of correlations and the relaxation time to the invariant density. 

2. The time scale when sticking to regular structures becomes important will be denoted by t*. 

3. In presence of external noise correlations are destroyed on the time scale t c ~ 1/c 2 , where a 2 is the variance of 
the noise. 

In order to observe the exponential decay of correlations it is required that 

tchaos <C t c , t* . (87) 

Destruction of sticking to regular structures requires t c < t* . In the regime 

tchaos "C t c < t* (88) 

the decay of correlations is expected to be similar to the one of hyperbolic systems and one should be able to explore 
it in the framework of the method that was outlined in the present review. Application of the method to specific 
examples is of great interest. 
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